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Abstract 



We prove some existence and uniqueness results and some qualitative properties 
for the solution of a system modelling the catalytic conversion in a cylinder. This 
model couples parabolic partial differential equations posed in a cylindrical domain 
and on its boundary. 

1 Introduction 

A gas containing N —1 different chemical species is flowing through a cylindrical passage 
with a parabolic speed profile. Chemical species are diffusing in the cylinder and are 
reacting only on the boundary of the cylinder. 

We investigate the existence, uniqueness and qualitative properties of a non- stationary 
system of partial differential equations describing the evolution of the concentrations of 
the — 1 chemical species and of the temperature, both with respect to the time variable 
and along the cylinder. This model of catalytic converters starts with the contribution of 
Ryan, Becke and Zygourakis [H]. However, we have added an axial diffusion term on the 
boundary in the present model, see [Tj for a more detailed description of this model. 

Due to its internal symetry, the cylinder may be reduced to the domain Q = [0, 1[ x 
]0, 1[, the boundary of which is S = {1} x ]0, 1[. The concentrations (resp. the tempera- 
ture) inside the cylinder Q are named Cif, i = 1, N — 1 (resp. Cnj)- The concentrations 
(resp. the temperature) on the boundary E are named Cis (resp. Cjqs)- The problem is 
written in a normalized way as: 



< 




(1) 



for i e {1, . . . , N} and with {z, t) = sup {Cis, 0). The initial and boundary conditions 
are: 



Cif{l,z,t) = 



dr 



{0,z,t) 



0, 



da. 



(2) 



(0,t) 



0. 



OZ OZ 

The functions Tj, i e {!,..., A^}, are supposed to be Lipschitz continuous: 



N 



2 I 
/is 



/i=l 



and verify the following hypotheses: 
(HI) Vx e (M+)^ : (xi, . . . , xjv) > 0. 

(H2) If at least one of the Xi, 1 < i < N , is equal to 0, then Fj (xi, . . . , 0, . . . , x^) — 0. 
(H3) For every {x,y) G {B+f x (M+)^: 



AT 



Si— (fj (xi, . . . , xjv) - (yi, . . . , y^)) {xi - Vi) > 0. 

lis 



i=l 



Remark 1.1 We observe that: 

1. Hypotheses {H2) and {H3) imply: 



N 



Vx e {R+f : - J2 ^i—'^i {xi,..., xn) X, > 0. 

i=i 

2. We have: 

Qf (1, z, t) = C,s {z, t) ; Qf (r, 0, t) = Qo (r) ; {z, 0) = (z) ■ 

In order to ensure the continuity of the concentrations in Q and on the boundary E 
at z = 0. we must have at t — 0, Cif (1, 0, 0) = Cj^ (0, 0), which implies : Cjo (1) = 
CisoiO) ;Q,(0,i) = Qo(l). 

2 Existence of the solution 

We establish the existence of the solution using the diagram: 

r 

Cif Cif 

Cif: 



Indeed, the proof of the existence is decomposed in two steps: 
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1. Existence of a solution in Q (given the solution on the boundary S) and existence 
of a solution on the boundary E (given the solution in Q); 

2. r = \I^o<|)isa contraction in some appropriate functional space. 



2.1 Preliminary results 

We have the following continuous embeddings (see [3', p. 103]): 



Because ^^(1, z,t) = j- —^"^ (1 ~ dr, we can write the problem as: 
dCif ^ I d f dCi 



r(l-r2) 

1 ri 9c,^ 



(1 -r^ 



dz 

da.. ^ d^a 



(3 l^lr'^f 
^ r dr \ dr 



lis rl 9Cif^ 1^ ^2 



r (1 - r2) dr + (5iri(C+ 



(3) 



with i G {1, . . . , A^}, and the initial or boundary conditions: 

(0,^,t) 



af{r,0,t) = Qo{r), 

■ Cis {z, t) 

■ 0, 



afii,z,t) 



dz 



dr 

Cisiz,0) 



0, 

CisO {z) 

0. 



2.2 Existence in the cylinder 

Assuming that the Cis are known on the boundary and performing a change of function 
in order to have homogeneous boundary conditions at r = 1, we obtain the following 
problem in which we omitted the time variable t: 



'l-r' 



duf_ 
dz 



n 1^ fdUf 

r dr \ dr 



Uf (r,0) 

Uf{l,z) 

{0,z) 



9uf 
dr 



dus 
dz 
uq (r) 

0, 
0, 



(4) 



with i G {1, . . . , A^} and 



Wf 



* (Ci/, . . . , Ctv/) 

Wf - Us, 



Us 



* {Cis, ■ ■ ■ , Cns) , 

diag (pif, ■ ■ ■ ,pNf) ■ 



We set 



«G(L2(0,l)f |^G(L2(0,1)) 



w 

Wo = \uem.{0,i)f 

and let Wq be the dual space of Wq 



du 

dr 
du 

dr 



N 



G(L2(0,l)f ,«(!)= 
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Definition 2.1 Assume that uq E {Ll (0, 1))^ and ^ {z) G (L^ (0, A function uj 



^N 



IS a 



weak solution of ^ if and only if Uf G L^(0, 1; Wq), ^ G L'^{0, 1; W^o) and 



or 



every (p G I^^(0, 1; Wq), we have: 
.1 /duf 



duf dip 
dr ' dr 



drdz 



Proposition 2.2 Let Uq and ^ be as in the precedent definition. Then there exists at 
least one weak solution of This solution verifies: 

duf 



'^°°(o,i;(LL ,,(0,1)) 



i-(l-r^) ^ 

2 



dr 



du 



(L2(0,l))^ 



dz 



< c, 

dz < C, 

< c, 



(5) 



0,1:W' 



where C is a positive constant which only depends on the data of the problem. 

Proof. We use a Galerkin approximation of u/ for which we establish the three above 
estimates, and pass to the hmit in order to prove the Proposition, see PP for the details. 



Proposition 2.3 Under the conditions given in the Definition \2.1\ the solution of ^ is 
such that: 



Proof. From the inclusion Wq C {Lj (0, 1))^ C Wq, and the fact that m/ G ^^(0, 1; Wq) 



dUf_ 

dz 



G L2(0, 1; H^o), we deduce the result using Proposition 23.23 of [Uf p. 422]. 



Remark 2.4 In order to take into account the time variable t, all the above expressions 
V {.) E H have to be understood as v {.,t) G (0, T; H). 



2.3 Existence on the boundary 



Assuming that is known, we have: 



dz 



— o 



dt dz"^ 
Us (^,0) 

^ dus , 



Us 

dz 



5r«) - j; (1 - r^) dr, 

0, 



(l,t) = 0, 



(6) 
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with: 



r 

6 



(ri, . . . ,r7vj 
■ ■ ■ ,6n) 



diag ( 



Is, - ■ ■ , Ons) ■ 



Let 



i/^(0,l) = |«G(L^(0,l))^||^G(L^(0,l))^} 
and let if* be the dual space of H^. 

Definition 2.5 Suppose that belongs to ^0, T; ^0, 1; (L^ (0, 1))^^^ and that Uso 
belongs to (L^ (0, 1))^- A function Ug is a weak solution of (0j if and only if Ug G 
L2(0,T;/7i(0,l)), G (0, T; if* (0, 1)), satisfies u,{z,0) = Uso{z), and if for all 
ipeL'^ (0,T;ifi (0,1)), we have: 



lo lo ( dzdt + jl (Os^, ^ ) dzdt 



dz dz 



lo lo (^r«), ^) dzdt - /g ( T/^, <^ ) r (1 - r2) rfrdzrft 



9m ; 



Proposition 2.6 Lei anc? u^o as in the preceding definition. Then there exists at least 
one weak solution Ug of (0). 

Proof. We prove the existence of a solution for the linearized weak formulation of the 
problem due to Theorem 2.2 of |T, p. 286] and then we use some fixed point argument. 



Proposition 2.7 Under the conditions given in the Defi,nition \2.f\ the solution of |^ is 
such that: 

Us G (0, T; (0, 1)) n C (o, T; (L^ (0, 1))' 



Proof. From the embeddings (0, 1) C (L^ (0, 1))^ C H* (0, 1) and the fact that 
Us G L2 (0, T; (0, 1)), ^ G (0, T; H* (0, 1)), we deduce the result using Proposition 
23.23 of p. 422]. ■ 



2.3.1 r = o $ is a contraction 



1 /2 

Consider the mapping T : Uj Us ^ Uf, and let /i = supj {'Jis/Pif) / infj 9is. 

Proposition 2.8 1. The mapping ^ : Uf ^ Ug is such that if we suppose fia"^ < 4, 
then: 





dUs 


ff 

lo Jo 


dz 



dzdt < 



4/i 



sup 



a"^ ^ia^)se[i),i]Jo Jo 



T pi 



Uf 



{r,s,t)r{l-r^)drdt. (7) 
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2. The mapping : Us Uj is such that: 

rT p1 



sup 

se[o,i] Jo Jo 



\Uff{r,s,t)r{l-r'^)drdt < - 







Jo Jo 


dz 



dzdt. 



Let : 



Wf = (0, T; (0, 1; Wo)) n (0, T; C (o, 1; (L^ (0, 1))^ 
W, = (0, T; (0, 1)) H C (o, T; {L' (0, 1))^ 



and 

Wg = {{u,v) eWf ^Wp\u (1, z,t)=v {z, t)} . 
Theorem 2.9 Under the hypothesis {H3) the problem (Qp admits a solution if: 

2 i \PifJ 

This solution belongs to Wg. 
Proof. Using ((7j) and ((HI), we obtain: 



(8) 



ITT II < ./ 



! Il^/llf ■= sup / / \\Uf\f {r, s,t)r{l — r'^)drdt. 
f se[o,i] Jo Jo 



The a which minimizes the Lipschitz constant is given by = 2//i, which leads to 



\\Uf\\j < ii^fe 



Ui 



/2. This proves that T : Uj ^ Uf is a. contraction if and only if 



3 Uniqueness of the solution 

Theorem 3.1 Assuming that the solution of (OJ) is smooth enough, the system (Op has at 
most one solution. 

Proof. We suppose the existence of two couples of solutions {Cjj.C}^ ^ and 

(C?f, ^£),=i,...,^ of (IID> and define: W,f = C}f - C}^ ; W,, = Cj, - Cl. We multiply the 
z-th equation of the system verified by Wif by rWif, integrate on [0, 1] x [0, 1] x [0, T] and 
use the equation verified by Wis on the boundary. Thanks to {H3), one can deduce that 
Wif is equal to zero in the cylinder because Wif is equal to zero at the inlet {z = 0) and 
at the outlet of the cylinder {z = 1) and its partial derivative with respect to r is equal 
to zero too, which implies, thanks to the equation verified by Wis on the boundary, that 
the partial derivative of Wif with respect to z is equal to zero in the cylinder. We also 
deduce that is equal to zero on the boundary. Wis is equal to zero at time T and at 
time because of the initial conditions. Because of T is arbitrarily chosen. Wis is equal 
to zero on the boundary S at any time. ■ 
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4 Qualitative properties of the solution 

4.1 Nonnegativity of the solution 

Proposition 4.1 For almost every {r,z,t) m ]0, l[x ]0, l[x ]0, T[, and fori G {1, . . . , A^}, 
we have : < Cij {r, z,t) , < Cis {z,t). 

Proof. This is obtained multiplying the equations of by the non-negative parts of 
Cif or Cis, respectively. ■ 

4.2 Upper and lower bounds of the concentrations 

Proposition 4.2 1. Let6i = —1. For almost every {r, z,t) in]0,l[x]0,l[x]0,T[, we 
have : < Cj/ (r, z, t) < AiQ ; < Cis {z, t) < AiQ, with: 

Aio = max sup Cio (r) , sup Ciso (z) 
\re[o,i] ^e[o,i] 

2. Let 5i = 1. For almost every {r,z,t) in ]0, 1[ x ]0,1[ x ]0,T[, we have: < 
Cif (r, z, t) , aio < Cis {z, t), with: 



= min inf Cio(r), inf CisQ^z] 

\re[0,l] ^£[0,1] 

3. Let 5i = l. We have for almost every (r, z, t) in ]0, 1] x ]0, 1[ x ]0, T[ ; Cif (r, z, t) < 
ttioe^* , Cis{z,t) < aioe^\ with: 

Uio = max sup Cio{r), sup Ciso{z) ) , A = sup /cj 

Yre[0,l] zG[0,l] 

ki being the Lipschitz constant of the function r^. 

4- Let 6i = 1. There exist two positive constants a and b such that we have for every 
I in ]0, 1[ and for every T > 0: 

J^{C^s)Hz,T)dz < aT + b. 

Proof. The verification of these qualitative properties of the solution is essentially 
obtained multiplying the equations of (Q) by the appropriate non-negative or non-positive 
parts of the corresponding test-functions. ■ 
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5 Numerical simulation for the reaction CO + O2 

CO2 

In the fluid, we use the following discretization method: 

• If z is different of N (chemical species), we directly solve (0)i using the method of 
finite differences or that based on finite elements. This requires that the equation 
with i = N has already been solved in order to put the appropriate values of mat. 

• If z is equal to N (temperature) we evaluate the coefficients at the step before. This 
possibly requires the use of some fixed point argument. 

On the boundary, we still use some finite differences method or some finite element 
method (see Pj for the details). 

Let us take the following initial and boundary conditions: 



in the cylinder: 



CO (r,0,t) 
O2 (r,0,t) 
CO2 (r,0,t) 
T(r,0,t) 



0.02 
0.05 



500, 



on the boundary: 



CO{z,0) 
02(^,0) 
CO2{z,0) 
T{z,0) 



0.02 
0.05 


490, 



We have the following graphs at t = 0.3s, 12s, 24s, 36s and 60s. We observe that 
the CO and O2 concentrations are decreasing and that the temperature and the CO2 
concentration are increasing: 

- Temperature 

502 
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CO 



C!,dI|S 




Oo 



0.0498 



0.0496 



0,0494 



Q.fl492. 



1.0488 
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- C02 




' ■ ^ 1. 

e 154 0.6 qm-: 1 

Remark 5.1 The reaction ends after 54s with the following values at the outlet of the 
cylinder = l).■ 
CO(l,54) = 0.017777, 
02(1,54) = 0.048912, 
C02(l,54) = 0.002355, 
r(l,54) = 500.873497. 
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